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ABSTRACT 

A program was written for the IBM 360 system of the 
Naval Postgraduate School Computer Facility to apply the 
methods of singular perturbation theory to earth-moon 
trajectories . 

Verification of results against previous work was 
carried out and it was found that agreement could be 
attained only for energy parameter values of 0.707 or less 
No solution for higher values could be found. 

The analysis of three-dimensional orbits was then con- 
ducted within this restricted range to show the merit of 
singular perturbation theory as an initial design tool. 
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I . INTRODUCTION 



With the advent of space travel and of complex satellite 
systems, a revival of celestial mechanics has appeared in 
modern science and engineering. In its present form, 
commonly referred to as orbital mechanics or astrodynamics , 
it combines the theories of astronomers and scientists such 
as Kepler and Newton with modern ideas and technology so 
that man may navigate his way through space. 

Orbital mechanics problems grow in complexity as the 
number of bodies under consideration grow. In fact, they 
become so complex that only the problem of two-bodies may 
be truly solved. Solutions to forms of the three-body 
problem are accomplished under certain special conditions, 
but usually numerical approximations exist for these higher 
order problems. Szebehely [Ref. 1] states, 

’’Not only , is the two-body problem solved — and 
the meaning of ’solution* may be different for 
astronomers, engineers, and mathematicians — but a 
general understanding exists regarding this dynamical 
system. The problem of three-bodies on the other 
hand is neither solved nor is the behavior of the 
dynamical system completely understood. 

In order to cope with the problem of space navigation, 
many methods to determine the motion of a vehicle in space 
have been attempted. The use of two-body solutions has been 
used in the method of patched conic sections. Many pertur- 
bation techniques have been applied for two, three, four, 
and more bodies. The most common method of calculating 
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these trajectories is probably the use of the digital 
computer. The computer is used to solve the initial-or 
final-value differential equations of motion in numerical 
form and is often time consuming and expensive. 

Due to the difficulty in finding solutions to trajectory 
problems 3 Lagerstrom and Kevorkian in References 2 and 3 
formulated a method in which they attempt to treat the 
restricted three-body problem with a technique similar to 
that used in certain fields of fluid mechanics. This is 
the method of solution of singular perturbations by inner 
and outer expansions and represents an attempt to obtain a 
uniformly valid analytic expression for trajectory motion 
in the restricted three-body problem. The work of Lager- 
strom and Kevorkian is aimed at giving a tool to more easily 
predict the effect of initial conditions on the motion of 
a vehicle as it becomes influenced by a third body in its 
field of space. This would allow the determination of 
’’ballpark" data for use in more detailed solution programs 
at a considerable savings of time and expense. 

In this paper 3 an attempt to use the singular pertur- 
bation method of Lagerstrom and Kevorkian to study effects 
of initial conditions on lunar trajectories is made. A 
brief discussion of Lagerstrom and Kevorkian’s theory is 
followed by some numerical applications of this theory and 
the results are discussed. 
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II. THEORY 



A. EXAMPLE INNER AND OUTER EXPANSION MATCHING PROBLEM 

Before formulating the singular perturbation solution to 
the restricted three-body problem it is useful to look at 
a simple example of matched inner and outer asymptotic 
expansions. The example is basically that of Prandtl and 
has been used by Van Dyke [*1], Ames [5] and others as a 
classical singular perturbation solution. 

If an equation of the form 

/iX^o.Fbc.x'-- x">V o 

is written and a regular expansion is attempted, difficulty 
is expected since the equation obtained when = 0 is of 
lower order than the original. 

Consider the following equation 

a) 



with initial conditions, 

? - o and $ ' 

The exact solution to this equation is 

f(x:^ - - e- v >*J /[\-e /M \ + t X (2) 

If a regular perturbation expansion is attempted one may 
write 
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f sf. 4 4 ^ + • • • 

i; - ■ - - 

r*c ->*»•«/*»; -■ • • • 

substituting into 1 
matching powers of implies 

t: - \> 

so that by integration 

t, = k* 4 C. (3) 

The boundary conditions cannot be satisfied unless b = 1. 
If the boundary condition at x = 0 is dropped 

= (j-M 4- b x 

since 

l = Uo + c — * C. = l— b 

therefore 

^oCx,/'-') = bx + Cl-b') (if) 

This approximates the solution quite well except in 
the region adjacent x = 0 (commonly referred to as the 
boundary layer), see Figure 1. 
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Ames states that if terms are lost or boundary con- 
ditions discarded in the outer solution they must be 
included in an inner expansion which must be developed 
(i.e., this implies matching in the overlap region). 

The asymptotic matching principle states that the m-term 
inner expansion of (the n-term outer expansion) = the n-term 
outer expansion of (the m-term inner expansion). See Figure 
2 for a diagram of the calculation sequence applying this 
principle . 

Returning to the present example, the first term is 

?o = U-M + 

for the outer expansion. The first term of the inner, 
expansion is found by a "coordinate stretching." The 
dependent and the independent variables may be stretched 
together or separately using physical insight or trial and 
error, however, the stretching must include essential 
elements omitted in the outer expansion procedures and must 
be matchable to the outer expansion. 

In the present example one may desire to stretch only 
the independent variable such that, 

fc*,>o = 

where 

X = Vsc ^ 
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Using this in the initial equation 



J 1 ? . S(/±l S’C^ b 

7x- ' r jF * jz- 



( 5 ) 



with 

F<°> - o , F Ols) =■ | 

where 



3(/*)o(X - <4 /t 

The highest derivative was lost in the outer problem 
2 2 

therefore d F/dX must be kept for the inner problem. This 
implies 

SC/A//+ ± CO as A - ° 

but 



^ ~ ° as ^jla. — * o 

implies that the solution satisfying the inner boundary con- 
dition is a multiple of X and cannot be matched. Therefore, 
the only possibility is that sCy) = a constant as 

* 0, and this constant is assumed to be 1. 

This gives 



dif 4 HL 
4x* 4X 



b^ 



( 6 ) 



with 

FCcO = O , FO/^ •= I > SLy>> - /A. 

In the inner variables X = x/y and F, the following expan- 
sions are produced: 

V - F o *t V % -4 r -U + ... 
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dl , 

AX 



- 

AX' 



43l 4 y, 4El + ^ A£. 4 • • • 
AX A AX ' AX 



£h_ 

AX" 






!£' 

l-v l ' 



Ax 



4lli + • 
Ax’ 



which when substituted into 5 yield 

+ — = © and lr (o') - o 

AX' AX 

and finally, 

Fo ex') = AO- e~ X ) (7) 

where the outer boundary condition has been discarded. 

To carry out the matching, 

1-term outer expansion = (.1- t % 
in inner variables tjA-X 

expanded for small jla = £l-t) + 

1 term inner expansion = O-k) 

1-term inner expansion = A C I - C. A ) 
in outer variables = A ( \ - 
expanded for small = A 
1 term outer expansion = A 

matching first term outer to first term inner gives 

A= l-b 

therefore 

To - Cl-b)Ct - e" x ) 

so that now 

P r V • + ^ + y»y Pi, + * * * 
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F = C l-t)C I - ) 4 Terras 

^ ^ + }^ X ^ ^ 

£ - Ci4)4-1d% + Mo MOM tbri^s 

to find the second term inner expansion 

- b 

MX* 

F, = C (e - * + 0 + tx 

therefore 

V - O-b)Ci-e" x ) + / ^L^(.e~ x +0 +t X]+ ••• (8) 

Since the exact solution for the outer variable is 
known, the matching from 1 to 2 may be skipped as it simply 
yields, 

I -t = I - b 

Going now from m = 2 and n = 2, 

2-term outer expansion = O-b} +■ b ¥■ 

rewritten in inner variables = £l-fe)+ 

expanded for small y 
2 term inner expansion 
2-term inner expansion 
in outer variables 
expanded for small y 
2 term outer expansion 
rewrite in inner variables 



= CI-W+ 

= C i-fc) + y- b X 
= 0-4)0 - e*) + / ^\c.(e y -i)+'=x] 
= 0 -t Xl - e 1 ^ ) 4 (gKi) 4 b t ^ 
= 0-t) -/-C4t.X + • •• 

= 6-4) + tx 

= O-b)^^ _ ^ c 
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matching these gives 

bx - — C 

C. - o 

therefore 

$ = (\-»=>) +■ b % 

f - Ci-tKi + (9) 

This form of solution is necessary when no single 
asymptotic expansion is uniformly valid throughout the field 
of interest in a problem. This problem may arise through 
factors such as two characteristic lengths .as those of the 
boundary layer thickness near the body compared to the chord 
length of the body. 

This long sample solution is presented above so that as 
the solution of the earth-moon problem is discussed an 
understanding of the procedures involved can be felt even 
though the actual solution is too lengthy to be presented 
in detail. The procedures are the same in the orbit problem 
as in the example just discussed but more complicated to 
carry out . 

B. NATURE OP THE PROBLEM 

With the method of inner and outer expansions previewed, 
attention is now turned to the trajectory problem in earth- 
moon space. It is beneficial to observe first the general 
problem which must be studied, the three-body problem. 

The three-body problem is one aspect of the n-body 
problem of mechanics. In earth-moon space the problem may 
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be depicted as in Figure 3* This figure represents the 
restricted three-body problem which is the second simplest 
of the many-body problems. The restricted three-body prob- 
lem, however, is still unsolvable in a closed analytic 
manner. The only n-body problem which has been solved is 
the two-body problem. It is therefore of interest to 
observe the nature of the three-body problem and determine 
what information one may find from determinable quantities. 

Using the system of three masses depicted in Figure 4 
it is possible to define three position vectors in space 
by r^, and r^. Now, if the conventional methods of 

vector mechanics are applied to the system and the necessary 
algebraic manipulations are carried out it is possible to 
determine ten of the necessary eighteen integrals of the 
total solution to the three-body problem. 

Writing the equations of motion with reference to 
Figure 4 





6 ] ^3 

_ 3 



( 10 ) 







+- Cl w\ 3 Yv\ x 



using the fact that 



and 
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and adding the three equations 



V* I f t t -t YV1 S f 



^3 - Ci Cy* * v* x r> * 



4 VW , VV1 , * T VvA^VVX 



- VA, vw, ^i 

A.\ 



W\ 3 W\ 






which yields 



0 



(ID 



which is the conservation of linear momentum. 

This vector equation integrated twice then gives six of the 
ten available integrals. 

Now using the relation 

0 * - ?t * A - ?! * ?» 

on equations 10 the following equations are obtained: 

kw. (?, } c c-i (j K xv- t ) + a y* (?, * a, ) 

^ ^>3 



KVU X.?; ) r^7 C^t. x? 3 ) 

Vx, fv J 

(7- 3 X.?-, ) +- Ci MtVAT' ^ 

fj, Oi 

If the vector relationship 

* f y - - fj X fc 

is applied and the equations are added the following equation 
is found: 

2 M : O'c * 'c ) - Ci £ 'jc?, x K) -La e - Vc \ 
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which reduces to 
noting that 

i t? xe ) = rx.r x7 = ext 

Cr+fc 

the equation 

21 vvv; dr, xtc) - o 
may be rewritten as 

J >*: [7; * Tt. ] J(t - wu L^c x r c ] 

which finally yields 

2 m; Hr; X XT; 1 = C.o».T (12) 

which is conservation of angular momentum. 

This allows the determination of three more of the necessary 
integrals for the three-body problem. 

The following series of algebraic equations will lead 
to the determination of the final integral for the three- 
body problem. Writing Newton’s law and operating on it 
one obtains: 

Ft - VA; 

F;. • ^11 - yvu • fc- ) 

Jtt 

Jf 4 

\ . if* 

Using this operation on Equation 10 
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z. 



w\, v, 2 



zz A 







\ q m 



^ - Ar v 

r x \ 



i- wixV'v 7 ’ = Cl kvuimJ 4 6\Vv\ z ^^\ r 2lJ^L 

J J ^ 



1 v*i I V - c, w 3 r*» f rii * j— I- 4- VAiWi, \ Y ?^ , cAr3 

J A? J AV 



the equations are now added to find 



- •£ tr- r 



J 



Af ^ 



-fC, WAl^ 




r — Ci wu 



i 




so that 



- ‘ Ci wi; y* ; 



2(4 yn - ^ — ^r*-— 1 ^ = co« 5 x ( 13 ) 

v> \ 

This is the tenth integral which represents the total 
conservation of energy. With eighteen integrals necessary 
for solution and only ten integrals that are defineable 
further pursuit of an analytic solution of the three-body 
problem is impossible. 

If a restricted three-body system is assumed one further 
integral may be found to offer new information. This 
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integral, called the Jacobi integral, does not solve the 
problem in analytic form however. 

The following two assumptions are adopted to formulate 
a restricted three-body problem: 

1) m^ « m^ and m^ 

2) m^ and m^ are in circular orbits about each other. 

With this preliminary introduction to the three-body 
and the restricted three-body problems concluded attention 
will now be turned to a singular perturbation solution to 
the restricted three-body problem formulated by Lagerstrom 
and Kevorkian [2] . The development will be outlined in 
notation of Reference 2, however for more detailed dis- 
cussion and mathematical details the reader is referred 
to Ref. 2. . 

C. FORMULATION OF THE RESTRICTED THREE-BODY PROBLEM 

With a brief discussion of the general three-body 
problem completed and with an example problem of matched 
inner and outer expansions discussed, attention is now 
turned toward a discussion of Lagerstrom and Kevorkian T s 
work in the application of this expansion solution to a 
restricted three-body problem in earth-moon space. 

The planar restricted three-body problem may be formu- 
lated considering the notation used in Figure 5. In this 
figure x-y is an earth centered frame of reference and x- 
is a moon centered frame of reference. 
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Motion in the x-y system may be considered to be 
Keplerian about the earth with only small perturbation 
terms due to the moon until the particle comes into the 
influence of the moon’s gravity field. At this time the 
perturbation from the moon is no longer small, hence the 
second characteristic length which results in the singular 
perturbation problem and allows the use of the boundary 
layer ’’stretching" type solutions. The effect of the moon 
destroys the uniform validity of the Keplerian solution 
about the earth. If the mass fraction of the earth-moon 
system is used as a scale factor one may discuss the motion 
near the moon in terms of the inner or x-y stretched 
coordinate system. 

Considering the outer problem first it is possible to 
write the planar equations of motion in x-y space directly 
as : 

+ ci -rO ■ = y- \ 

cH*- r’ 

( 1*0 

Ui-^1 - 

where 









- V 



and 










If the terms on the right hand side of the previous 
equations are set to zero, one has the Keplerian or two-body 
equations near the earth. The Keplerian integrals will 
change little near the earth. These integrals are those 
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of energy , angular momentum, and position. Recalling the 
equations for the three-body problem discussed earlier one 
has , 

2 < 'Oi ) - Co T 

and 

2 (.’t ^ !2l j ) ” Co MiT 

Rewriting these in two-dimensional coordinates in the 
notation of Lagerstrom one has 

t[ + c \ i ] - -Co-XT = U e 05) 

and 

XV^ - ^ x = C_o«*7 =~ (16) 

Two equations specifying the position of the orbit in 
the x-y plane may also be written as 

'fc - X € * + Cl — r 

(17) 

V = 4- CV-^^/r 

These are written such that if p g = 0 the motion is 

symmetric about the x-axis. Also if q e > 0 the x-coordinate 

of the perigee is negative for J- 0 and always non-negative 

for £ = 0 . 

e 

If the initial conditions are specified in terms of 
these integrals, (which will be done shortly) the motion of 
the particle in the Keplerian orbit may be determined. 

The launch is considered to be at perigee, hence 
P e = 0 is given. 
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The angular momentum is defined also to be rv cos <f). 

If perigee Is considered <J> = 0 and cos <f> = 1, therefore, 

& e = rv. 'If this Is non-dimensionalized with respect to 
the angular momentum of the earth-moon system (moon’s 
period) the non-dimensional angular momentum may be defined 
as , 

v rjs 

A " + 08 ) 

The non-dimensional energy is defined in terms of a 
non-dimensional potential energy and a non-dimensional 
kinetic energy as, 

f) 1 = + (19) 

Since h , J^ e , p e determine q e the particle initial con- 
ditions are specified in terms of h , P e ,£ e > and q g . 
Writing these in their non-dimensional form, 

k = -/»'> , ft-. 04 e?o and i - oCy?'"> 

one Is now able to specify the initial conditions of 
velocity, altitude, and position in terms of initial 
Keplerian Integrals. The following Keplerian relationships 
at perigee are known: 

— x P = r P = aL\ -O 
>< = o 

V = £-)*■»( V - t') 

p 3*- 

Replacing h by -p , i by u 2 A and p by 0 at x = 0 and 
e e 

carrying out some algebraic manipulations the initial 



values for the coordinates and velocities become: 

^ ! i _ _ . __ 



l Cl ) 



^ = ° 



qL>i - O 



( 20 ) 



At |-/^L 0~/*) ' 1 



where 



and 



A t 

a - g-/*> /"z p z 
e= \ \-l*-}*p'>?kt-r?]\' 



It is also beneficial to transform the motion equations 
to the independent variable x rather than t. One trans- 
formation is shown below and the rest follow similar steps 
to produce the transformed equations. Writing 

- K - 



At 1 ck* 



A Ax _ j. \ A 'MAk'|_L 

- jfc J5 - Ax lA+Ux I t' 



finally 



- A l tr] t' 

■ 'W F 



= V -c' 

At x (tV 



Vt" 

(t')' 



2 2 

A similar development may be used to transform d x/dt , h 0 , 
and so that these equations are finally available in the 
form: 
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-L 4 0-^0 ^ = y- $ 

£ 

f, - f?' ^ 

jk. _ Al'ii' - t^L^LHWl 

37 r \ r 



(.2 0 



Jt* 

$ ' il ' u ‘‘^" 

D. OUTER EXPANSION 

Attempting a solution of the equations of motion for the 
particle in earth-moon space in a power series in one has; 

+ ■■ • ( 22 ) 

ttx,^ - -to l-Xjh^ + -V - - • (23) 

The terms 2 y^ and t Q are the exact solutions for the 
equations of motion with the perturbation terms on the right 
hand side set to zero. In common with celestial mechanics 
notation, these terms are written as. 



where 



-to 


= G 3/l [tU)-£SwE O^] 


(2*1) 




'V^ -Ur-oEL*) 


(25) 


EC*) 


= -Ce-^ a ] ,/a 










a - 






e - 


[i-zyp^^/ci-^y z J /7 ~ 
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C2(S) 



If and are further expanded as 

toCyjyl — too £*) 4- c>o 4 • • • 

^(-V^ - Vj-*) + l-M^U • * • (27) 

where by appropriate expansions in powers of y< one finds. 



to.Cx) o 1 


[$| w-'zpCx- Cx -p'-x'V'’'] / 'Z- 3/l y0 3 


(28) 


to , co - 




(29) 






(30) 


4/*, <-*) 


= -XpCr-x ■*. >0 )( \-zp^)jzp LzjCx 


(3D 



It is also necessary to insure that the moon and the 
particle pass within the scale distance y< of each other. 
The phase constant expansion which includes this condition 
in the analysis is, 

T=T. + / /x T,i, + /-T, (32) 

where 

To - too CO 
~l1 - - (-■'I 

and is an arbitrary term which assures that the moon may 
be allowed to interact with the particle as desired. 

The perturbation terms t^x) and y^x) may also be 
determined, the details of which are long and tedious and 
may be found on pages 689-690 of Shi and Eckstein [6] . The 
final form of these equations is : 
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( 33 ) 



where 



and 



f * 

£, (*) - \ 4° ( -'0 L - **<-*V\ A=t 

U 1 ( - x ^ c * | u t -too ~ (3*0 

«c-?) = -[e" t 6-/J l 4) , 'Vz ( 35 ) 



^HluTfUrp] 



(36) 



also f (z) and g (z) are given in terms of the moon's posi- 
o o 

tion coordinates as, 

£os tB') — \ ol — ^ / r —i 

Tv'/ ~ -T.l 

' 1 1 + 1 - z. t coi -t<ti \ 1 



(37) 



-- -T.I-=g- _ ^ _ T ^ ( 38) 

| I -t -fc 1 -- 2-i. CoS Lt-x-fi) -To] ^ V 

If t^(x) and y^(x) are now expanded asymptotically as 
x -*• 1 the following expansions are produced: 



n- 3 4 



■£,60- £ Z-Zf 1 -) ^OyO -x) — YCp') + OUl 



(39) 



<2.1.0 = 4- S(fO 4 00) 



where 



Y ^il' U*>L°<co-^) + < 577 ^^ 4- 



Q 



('JO) 



(^1) 
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(H2) 



with 



Q - O + 



vV-L 



(*3) 



and 



U - I / t do (-0 



(M) 



With the terms t , t, , y n , and y 2 y now available the 
O 1 -L 

outer expansion for t and y is known and is ready to be 



matched near the moon. 



E. INNER EXPANSION 

As the vehicle approaches the moon it comes into a hyper- 
bolic encounter with this body. The hyperbolic encounter may 
be expressed in terms of inner variables. These inner 
variables are written as, 

X = x- io>6fc -T ) - _ y <,i.o6t-T ) 2 - i -Ti -y- ? 

) \ ^ J 

where y. is the scale factor for the inner variables and t 
is a constant which must be determined. 

The inner variable equations of motion are written as: 

J z x + £ = _ ±c^2 ^ +uO-rt <;o ^ (^5) 

Ai' r 3 i'e r r 



4l5 + 



S 1H± (1)6) 

V'e ' 
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where 



Vc = ) -4- X T f c\ t>/ -4- y}- T 1 

y CZ- ?r) + 1Z -T 

In the hyperbolic encounter one may express the Keplerian 
solution in parametric form as, 

/ - a(.e -Co*WuO^0 T s/~£ (^7) 

4 - ate-cosku^s/-^ t cc^B- (^8) 

^ = 5. % [cs/Aia — ^ (^9) 

The Keplerian integrals previously discussed expressed 
for inner variables are 




% = ~1 A £ +- */r 



These integrals are related to the hyperbolic elements by, 

a. = '/^ , e. = /h I^T 7 ” 
f = -e , % = -e. c-°=>S 

To carry out the matching of the inner and outer expan- 
sions the values of y and t as functions of x along the 
approach asymptote must be determined. These functional 
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values may be determined by allowing u As u » the 

following functional forms of y and t are determined to be: 



s- at* - + 



(5D 



1 _ x —k • r 2/r ) ’ z ~* 

- -QT 4 n u , e 



+ o C9 ■') 



(52) 



The U 1 and terms are the velocity components on the 
approach asymptote. 

, i - QJLA 



V, - ~ ^? S I 

l-f TTC^ X ^ 



where 



UP + v/p- = zTR 



and A and B. are, 

/U - 



and 



z - -p 



If the negative y-intercept is considered to be K^, 
that is to say y = - for x = 0, 

K. ■= M -Ky , ( 53 ) 

o, 



and 



A - K, U, 



(54) 



U 11 and V^, the velocity components along the departure 
hyperbola, may be written also 

U, ^ C-Z-V/H 

|+ ZX. 1 

v/ - - y 1 +• 

" ' l+z.'KT.' 1 ' 
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F. MATCHING OF THE INNER AND OUTER SOLUTIONS 



The matching of the inner and outer expansions must be 
carried out at this time. Returning to the example dis- 
cussed, the principle behind this matching may be recalled. 
The only difference in the two problems is that in the 
present case, since only the first term of the moon inner 
variable expansion is needed, it is only necessary to 
rewrite the outer expansion in inner variables and then 
carry out the matching. 

For the particle near the moon the x-coordinate is £ . 
This occurs when t = T . If the x-coordinate is specified 
to be 




this implies that 




The phase relation was written as. 



T= To t ja w *-Ti U * y-T, - To-t ■v 



( 55 ) 



but with t = Tq the equation may be written as 



'i- / =■ * { lz. 



so that 



“X - V- -V oco 



with 




X-l 

r 
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therefore 



7 * *■ 



r 



If the outer expansion is written as 

t " £ooOO+- yuHt oi <>) +-fciOc^ + • * • 



( 56 ) 



^ c y.' 1 *- ^nJ^>o) +• A- • - • (57) 

and replacing 

TC “ 1 +-^- (-X - TliXj-j. - ) 
t - zfo»c<> +y- [l v -Ik? ')U' 1 +“toi CO + (| + v/' _ ') J,1_ Ao^ 

+ Y(p)] +- °ty-} 



too CO=To A-* B - t - ~To — y- t 

r 

- C>< +■ to, o) +■ Ci-hjO W + 

+ Yip') - T + OL 0 



also it is recalled that 52 was written 

- * 4 - 1 - — 

-ur • /o r<l^e 



c _ a . - 3/ m - 2.X +- o ( v -• ") 

-6= a 






The term (1 + U 2 )~^ / ' 2 l°^ y- does not occur in the later equa- 
tion. It must therefore be cancelled by some other term. 
Since the only arbitrary term in the inner expansion for t 
is 5 a term in must cancel the log jjl term, hence, 

't =- 0 + U l -y 3/ ' L E, ( 59 ) 



and 

t =U~'* -f-O+u*") ’’A.oyl-x) 4Y(p)-2; +oc0 ( 60 ) 



3 1 * 



y may likewise be determined to be, 

n - ib- _ -to i O') +Ti -£ tot-0 

\ ZU u 

Comparison of 51 to 61 and 52 to 60 shows 

o,- u 

y, 1 

If the equations are written such that, 

" Yf K - K, + ocf) 

1 U, 

and 

C\ = — K - S- + oCiN 

where 

-K - Tli --to.CO-f-Ti + S(^> -Y(p') 
tu 



H> - L >K - 



Then it is evident that 

K. - K 

Rewriting t as 

£ + |oyd-x) + L + oO - ) 

where 

L- toy Li) - Ili- + YCf} - C", 

and 



I - ^ -4 

0, 



- Vz. 

(X 



U 









-f ° Cx" 1 ) 



then 



( 61 ) 



( 62 ) 



(63) 



(64) 

(65) 



(66) 



(67) 
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L~ -A- i- a. 

O, 



FT 



\c 



°( Ut £ ' 



which completes the equation matching process and determines 
the unknown constants. 

These values are used in equations 50 to obtain from 
substitution 



2^ = HU* 1 " 

i - ku 

^ = Ku , -fCi + i;y lk ( 68 ) 

| = ku - v c i + 



The term may also be completely determined to be 

% - 6..0)- Th- -hYCf) + L - a’^ocj. - (69) 

The total or composite expansion solution is the outer 
expansion plus the inner expansion minus the common part. 

It may be written as; 

^ +■ + ^1^65) - ol* C* ) ] + °<-y) ( 70 ) 



t - toCxj^ + j j>i.|(-x'> + y [jE £*) - + ocy) (7D 

where 

oL* = - 0, *.~§ 'i/Ui (72 ) 

£>* - (>^ - i e Co^)/lJ, 

.The solution to first order after moon passage may be 
obtained by a similar matching procedure as that used above. 

Before leaving the analytic theory it is wise to note 
that this procedure of inner and outer expansions may be 
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carried out to higher orders of p. which increases the 
accuracy as the value of the parameter u increases. 
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III. THE COMPUTER PROGRAM 



The computer program used In this paper was written for 
the IBM 360 system at the Naval Postgraduate School Computer 
Facility. It has been run successfully on the OS/CMS time 
share remote terminal and the standard batch processor unit. 

The program is a straight forward computation with the 
use of Scientific Subroutine Package programs from the sub- 
program library. These programs are used to evaluate the 
integrals which appear in the solution near the moon. 

All computation is carried out in double-precision 
notation to maintain as much accuracy as possible. The 
integration routine is a double-precision sixteen point 
Gaussian quadrature from the SSP library in the Fortran IV 
language . 

The necessary inputs to the program are the non- 
dimensional boundary conditions discussed in this paper. 

From these parameters the program calculates the necessary 
constants, then the integrals of the motion equations are 
computed for x = 1 (near the moon). With these values the 
program evaluates the motion in the vicinity of the moon 
and prints the necessary parameters. Control is then 
returned to the input mode and a new series of calculations 
may be carried out. If p becomes too large the machine 
proceeds to stop and ends processing. 



38 



The necessary equations, a test case, and the computer 
program are contained in the Appendix. A flow chart may 
be found in Figures 17 and 18 . All equations are explained 
in References 7 * 8 , and 9 and the same notation is used as 
much as possible. 
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IV . DISCUSSION OF RESULTS 



A. VERIFICATION OF COMPUTER RESULTS AGAINST THOSE OF 

REFERENCE 7 

A computer program to solve the restricted earth-moon 
three-body problem was written. The first test of this 
program was against the data presented in Reference 7. 

Figures 6, 7, and 8 show the agreement of the present 
results with those of Ref. 7. The results for y and <5 also 
agree with Table 1 of Reference 7* As a final check 
examples one through nine of Reference 7 were tested and 
comparison of results was good. 

No data of any form for values of p greater than 0.707 
could be made to agree with the published data of Reference 
7. The inability of the program to predict accurate results 
for this case of p greater than 0.707 was unresolved. The 
probable cause of the instability of the solution is due 
to the denominator of the f^Cx) term rapidly vanishing near 
this point. The equation for fg(x) was rewritten in 
several forms, all to no avail. The equation was also 
expanded in a Taylor's series about 0. This expansion gave 
some improvement in the prediction of y and 6, but the 
improvement was not sufficient to be used for numerical 
work. 

B. PARAMETRIC STUDY IN THREE DIMENSIONAL EARTH-MOON 

TRAJECTORIES 

The program used to predict the results of Reference 7 
was modified to include three dimensional orbit consideration 



for values of i = 0(1). The modification still led to the 
same results for planar orbits and a study for inclined 
orbits was carried out. Cases were run to predict the 
effect of variations of energy, angular momentum, position 
of the moon and orbital inclination on the space vehicle 
as it approached the moon’s position. 

The present program predicts all of the important 
parameters near the moon. The effect on changes in initial 
conditions can be observed in the changes in these near 
moon parameters. One of the most important of these 
parameters and one that requires knowledge of the other 
parameters is the parameter of closest approach to the 
moon. This parameter has been chosen to show the effects 
of varying initial conditions on the moon encounter. 

Figures 9 and 10 show how the initial launch conditions 
are related to altitude and velocity at launch from perigee 
The altitude limit on perigee shown in Figure 9 indicates 
values of lambda must be less than 3.0. Figure 10 shows 
the velocity at perigee. Since conventional launch vehicle 
cannot attain values of velocity much greater than around 
*40,000 kms/hr an outside limit on velocity was taken to be 
50,000 kms/hr. This, then implies that lambda should be 
greater than 1.0. For a given launch condition the values 
of rho and lambda can be related to the velocity and alti- 
tude at perigee. The perigee launch gives the angle of 
perigee to be 7T and time of launch to be t = 0. Hence, 
this defines the launch conditions completely. 
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Figure 11 shows the effect of changes in the initial 
energy, rho, relative the earth on the distance of closest 
approach. It is apparent that as the energy is reduced 
(rho becomes greater) the distance of closest approach to 
the moon is increased. The effect of orbit inclination on 
the distance of closest approach is also evident. This 
Figure predicts an increase in the distance as the inclina- 
tion angle of the orbit is increased. 

Figure 12 depicts the effects of changes in the angular 
momentum, lambda, relative to the earth on the approach 
distance to the moon. As the angular momentum is increased 
(this corresponds to an increase in velocity and/or perigee 
height) the distance of closest approach to the moon is 
increased. The effect of the inclination of the orbit shows 
up in this Figure as before to indicate greater approach 
distances for higher values of inclination. 

Figure 13 can be used to predict two effects which occur 
due to parameter changes. The effect of T-^ on the closest 
approach distance is due to the variation of the moon’s 
position at launch of the vehicle. For positive values of 
T 1 the moon and the vehicle interact very closely, whereas 
for negative values of T 1 the moon passes far from the 
vehicle . 

The effect of a shift on the distance of closest 
approach may also be seen in Figure 13. Smaller values of 
the apse angle allow for closer interaction with the moon. 



The value for planar orbits is independent of the apse angle 
and therefore remains the same. 

Figures 16 , 17 , and 18 show the effect of these same 
parameters on the distance of closest approach for orbits 
of only slight inclinations. As can be seen from these 
figures the effects of changes in rho, lambda, and are 
the same as before. As would be expected changes in the 
apse angle and small changes in the inclination angle from 
the planar case cause no significant change in the distance 
of closest approach as predicted in the planar case. 

Figures 14 and 15 show the three-dimensional geometry 
of the earth orbit and the moon passage. The inclination 
angle, i, discussed is the angle between the earth-moon 
plane and the orbit plane. The apse angle, to, is the angle 
from the node line, fi, to the apse axis. Similar parameters 
are shown for the moon passage orbit. 
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v . CONCLUSIONS AND RECOMMENDATIONS 



The method of uniformly valid asymptotic expansions 
appears to have merit as an initial analytic tool. The 
ability of this solution to predict values of parameters 
of earth-moon trajectories with such a small expenditure of 
computer time well justifies use of these expansion solu- 
tions. A quick parametric variational study could predict 
"ball-park" initial conditions which would greatly reduce 
the time of total solution of a trajectory problem. 

Observing the results of the analysis of this paper it 
appears that one might conclude that the closest approach 
possible occurs when the launch energy is as high as possible. 
The angular momentum relative to the earth must be low also 
for this closest approach to occur. The moon f s position 
which is fixed by T 1 and w must be such that it is just 
ahead of the particle as they come near each other. 

The only other parameter which must be considered is the 
launch inclination. The minimum value of closest approach 
is indicated to occur for planar trajectories. 

From this discussion one then concludes that minimum 
distance of closest approach occurs for a planar launch at 
the highest value of energy and lowest value of angular 
momentum possible. These are attained for near earth 
launch with high perigee velocities. 



It is recommended that further study should be conducted 
in the vicinity of the minimum energy solution region. This 
is the region of the previous mentioned instability and is 
very important in the study of low energy transfers where 
most of todays travel takes place. 

It would also be interesting to release some of the 
restrictions on the analysis in order to apply the solution 
to interplanetary problems. 

Kevorkian [8] has extended the analysis to higher order 
terms and a comparison if this result with those higher 
order solutions would also be of interest. 

Additional work is necessary before the present analysis 
can be applied to a full spectrum of trajectory problems. 

The present work however indicates that useful information 
may be gained from this form of solution for the restricted 
three-body problem and that this solution has merit as a 
potential analytical tool for initial design studies. 




FIGURE I. AMES, VAN DYKE BOUNDARY LAYER PROBLEM. 
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Expansion Expansion 




FIGURE 2. FLOW CHART OF MATCHING PRINCIPLE. 
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FIGURE 3. EXAMPLE EARTH- MOON HYPERBOLIC ENCOUNTER 
TRAJECTORY. 
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FIGURE 4, THREE -BODY GEOMETRY 
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FIGURE 5 , EARTH -MOON COORDINATE SYSTEMS FOR THE 
RESTRICTED 3- BODY PROBLEM. 
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FIGURE 6. ANGULAR MOMENTUM CORRECTION. 




FIGURE 7. ANGULAR MOMENTUM RELATIVE TO THE MOON. 
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FIGURE 8. DISTANCES OF CLOSEST ENCOUNTER WITH MOON. 
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FIGURE 9. PERIGEE HEIGHT VS. LAMBDA FOR VALUES 
OF RHO. 




FIGURE 10. PERIGEE VELOCITY VS. LAMBDA FOR VALUES 
OF RHO. 
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FIGURE II. INCLINATION VS. CLOSEST APPROACH FOR 

VARYING ENERGY, ALL OTHER PARAMETERS 
CONSTANT. 
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FIGURE 12 . INCLINATION VS. CLOSEST APPROACH FOR 
VARYING ANGULAR MOMENTUM, ALL OTHER 
PARAMETERS CONSTANT. 
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I ~deg. 

RGURE 13. INCLINATION VS. CLOSEST APPROACH FOR 
VARYING . PHASE CONSTANTS, ) ALL OTHER 
PARAMETERS CONSTANT. 
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FIGURE 14, EARTH ORBIT 3-D GEOMETRY 
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FI&1RE 15. MOON PASSAGE 3-D GEOMETRY 
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FIGURE 16. APSE ROTATION ANGLE VS. CLOSEST APPROACH 
FOR 3 VALUES OF i WITH p AND T, CONSTANT 
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FIGURE 17. APSE ROTATION ANGLE VS. CLOSEST 
APPROACH FOR 3 VALUES OF i WITH 
XAND T, FIXED 
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FIGURE 18. APSE ROTATION ANGLE VS . CLOSEST 
APPROACH FOR 3 VALDES OF i WITH 
A AND X CONSTANT. 
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FIGURE 19. PROGRAM FLOW CHART PART A 
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FIGURE 20 PROGRAM FLOW CHART PART B 



THE INITIAL PARAMETERS ARE 
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COMPUTER PROGRAM TEST CASE 



LIST OF EQUATIONS USED IN THE IBM SOLUTION 



The flowchart of Figures 17 and 18 shows how the logic 
of the computer solution is carried out. The following 
equations are those used in the solution on the computer. 



INPUT 



f = i * 



j 1 



'lx. 



ypyp </ t 

j y- DCnC ) \ 



y- 



r. o . o i z i jT £ 



yvi VK 

7*1 r ?kcs< 



(nr m * lO"'° K^/kCi- t>Av‘ 1 ' 

P&rzc zil J — v>avj 

X - Ofi-ftlT |MC »fi / AJ C-vLb 

<4p s £? Zyc«r/o/^ 



'EVALUATION OF TERMS FOR X = 1 

Tool = £oo(J) - [5r~- , -Zf’0~/P*)'*--z r Ci-^yU/z-’V 3 

TooPI T too O') = C 1 / 

Y'*«" = Cl- p i V /l - 

CAPTO = To = too Ci') 



CAPT.t-B ^ AUCou 

CAPT •= T =■ To+F‘"“.^+-^-— 

-- - [o -f‘1' 1 / i + i/o v «y'*._ x 

u = U = l/ too CO 

<3 = 3 - C> + U i -') V ' 1 - 
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EVALUATION OF TERMS IN x-y (MOON) SPACE 
H'&AE. = - U4vJ v ) 1-2. 

XLSMl = J =. +^(.I+0 7 -Si^ 1 l)/ 2. +OCV-S') -utrorTi - )-^ U,o)|/coic 

k £ E - | / 2JK 

PlB.AiZ-5 f « O I + (.ilO- 1 '"’’ 

<3^ 5 5 . jf _ OCl^r 11 ' 1 - 

EPSBAH.H e = £/ + z.fc 

~TAi/ saij-s 2r = 6. Cy-*- 6, t e *£) + ro-avuu £«*s i 4c ,,l + lit i 

7«<=8Ae-=& = t*~~' C~7 Cl+U l 'i"' L ) 
lSA.fi- ~l ^ S;„-' (- £ C/ + L> l )"VJ 
<3^b =-a = J**-* C-i/u) 

EVALUATION OF VELOCITY COMPONENTS 

U | * U. =• u 

VI s V, 

W I r \aJ | CL- O 

l)il = On : 

vu = >/„ =(!/e^iW'^ ,+ ^ ir -'1 S/ ^^ d+u "' )c ^ 5:Cc ‘^ 

Wit s w lv =r [V { i yi+^) s/ajX 
EVALUATION OF EARTH ORBIT PROPERTIES AFTER ENCOUNTER 

P5J = V - 4 »~ C/Vm+O/Uu] 

HEU = Jte„ - 1i 4^ii - Vt 
£>EE = ^e„ - -(1 /m + iVm + I 
?&IT £ f>e„ - £ l/u 
X(_etTB J| eii _ v/(i + I 
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CLOSEST APPROACH 

dr - 0/zW)& 4 '* zK -i ] 

COMPUTE PARAMETER CHECKS OF LAGERSTROM AND KEVORKIAN 
T>AP-A> - UCr-S) = UCCFIl - CC-xii) 



7 A#AZ = J + U 



T, - XV' 



EQUATIONS USED IN THE SUBROUTINES TO OBTAIN THE MOTION 
INTEGRALS / i« A 

T oo% 3 -2^r^ 3Ky 

Too?i s -took.) =• — p 7 -^)V'~ 

licpHi i*. = -[t' v 6 -pV/V 2 - + i i/z-^i 



f 



Fo? B U) = Co^-bo^ -Tol - * _ COS &,.<*> -Tol 

£f-f l 2 - Z.^Co^U: OoCt) — T O J S 

C -.01 = %U) r s '~ 2 _ -To 7 

£ f+ CoS C’t’ooCi') -To] $ *- 



CF “ = Vrr L + i 

cei " E 6 = f # ‘ \ i *j.te) [ pco + i 
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THIS PROGRAM USES THE METHOD OF MATCHED ASYMPTOTIC EXPANSIONS TO SOLVE THE 
SINGULAR PERTURBATION PROBLEM IN EARTH-MOON SPACE* THE ANALYSIS IS BASED 
ON WORK BY LAGERSTROM WHICH HAS APPEARED IN THE A I A A JOURNALS* 



oooooooooo o 

—'(MO-'vt IT'D oca O — CM 
000000 OOr~» — — 

oooooooooo o 
OOCCCOCCCO o 
cpccocccocrccccecco cn 
cl cl cl cl cl cl a cl a cl cl 
OCCOCOCCGC c 



ooooc ooooooooo oo 
o— oico>t *T(>Nrvooc'o — rr-in 

!f\)rorr r^r' vt <► cf vj- 

oococ ooooooooo co 
ooococoocccoco oo 
cocrcocccccccc'cccc'cc ccccccrc cc cc 
clclolclclclclclclclclclclol cl cl 
DCOCCQCCCCGOCC OO 



2 

3 



a: 

o 



S’ 

o 

o 

X 

CL 

< 



>- 

o 

o 

CD 



< 

o 



— o 

— C\J 

rvj>- 

? ► 
X — 

+ o 

HO 

si- 

X ► 

*- CL 

oo 

— O 



0 

1 
CL. 

* 

o 

o 

0 

CM — 

1 —— 

— o— 

— CfM 

H — O* 

— OJCM* 
I- — * VO 
CL —-ft Ol 
<1 — 00 CL 
Cl <\jX e — 

* CL CO I 

* w^-O 

I o 
II XO* c 
XorGCr 



o 


K 










_ - 


— oG — 


k- 




III 


UJ 


H 






COO 


OC 1 H eh 


CL 


21 


O 


O 


K 






#■ ► 


OO — CM CL 


< 


O 


CL 


X 


► 






* H 


U.CI — o 


O 


— 5 


< 


< 


O 






O <f 


CCL-K-O0 


* — 


1 — UJ 


K 


»- 


x 






— X 


00— 0—0 


— O 


o_j 




OO 


X”) 






— CL 


21 — OO rc, * 


OO 


sec 


«*O0 


— 


•* ► 






K-J 


CL« O* CC 


X e 


o 


►00 


o 


oco<r 






OCL< 


UJO-X * X 


— XfM 


ID cL 


NJ< 




►so 






CT O ► 


KCCO< 




H-CL 


1 x 


IZ 


(\l<3 UJ 






OOOO 


»IIJ 


XKO 


< 


o 


IDO 


zjf 






Ck- 


00 CM CL CL X 


— CLO 


_JX 


► CVOOI— — 


X X O 






C* CL 


»--* — * 


ooo e 


30 


XXUJoOk- 


► -oo 


o 




KO< 


X vo— — 


Coo co 


OO 


1 Xl/)X < 


HO ► 


O' 


oo 


oo 


<rocvo 


oc — 


-JCC 


<r ►oooocl 


xx— 


o 


CL 


OsO - 


k-O « — o 


o * 


c 1 


— 


UJ 


X CL X 




UJOOhO 


oo o (V— c 


■i + ■* 


(Jfc 


<CD25.U.d 


W\ 


o 




I— vtl 


X — — CM 


3 — . 




»-* UJ 


CUJ 


—IfOCM 


k- 


UJ | 


Of— *CL 


C— X — — 


— k O CM 


OO 


C-JCLLL 


oo 


<r<r<r — 


— 


X LL'O oVHCiUh— Mh- 


CL CCk-*- 


KUJ 


D«CQ 


I-KI- — 


co c\j O 


<C 


! ©fOU-k-k- 


OcOO-k CL 


oxcl* 




ID CL 


Ol- 


<«r«<0 


CLO— -k DZZJC'CLr O 


CC<3 



so 
<r — 

CL CL 

Oh 
O oo 

or. Uj 
CL CL 



ICL 0 — 01 

< oo— clz: 

— I — OO I— LU < 

k — «aj 
-051 CL CL 
X — I 

— -l—o CL I 

clx xiu 
xxxclo 



CCOhhl 

VXV.<«TCO 
hhhOCC o 
COLOOO •- •- -UJ 
*— « • — •— LO to to ^ 

LL!UJUjOCjC-> 

< << LLIUJLUIX. 
XXXcLCLCL — 

CO 



<r O • OCtOoO_JcO<OUO — I— -JO + 
CLOhOO <TOOICOVX o 
CmO ©<_J_J — CL OO II II O 

LDV ofMV<<UJ II — II K-OtM c 

hvOO ZZZh II I II © — — 

GU'IIICO'OCGH 0 — 0 — I- Cl — 
Cl « CL CL X UJ UJ CL CL — C Oh- CL < 
S't-ujri-i-rco ©CVCL ii < o ii 
O|IU.0-CXXOOCHH< o 

uc w oujaiui-h'->-uD o 
< 

CO 

o 

O 



ooooooooooooooooo ooooo 



67 



oocccoocooc oooooooo 

>nx ooOrrJr rvJ ir *c cr O r-i r\i<nn x- cr 

st <} iruniOL/M/Mnin in vO >o *o >0 >o *c O 

ooooococccc ocoococo 

OOOOOGOOOOC OOOOOOOO 

CC CCCDCOCC CCCC'CCcCCC’C: cC CD CO CC CC CD CC CD 

cLcLcLaacLCLCLcr'CLa' clclotclclclclcl 
CCCOCCCCOCC CGGOOOOD 
* 

■W vD 



<M 



oo OOOOOOOO 
>oxcr OHf v :f^st m 
t — T r~* X- 00 CC CG CO CC CO 

ocooooc ooo 
OOOC OOCoCO 
CDCCOL.‘CCCCCDCCCr.CCiCC 

<X Cf-O^C£CX CC CC CL CL CC 

COCOCCCCOL 



o 


ooo 


o 


vO 


Oi— IO 


>* 


CO 


aoo 


O' 


o 


ooo 


o 


O 


ooo 


o 


CC 


coco CO 


co 


a: 


CLCLCL 


CL 


o 


CCG 


□ 



o 


* 


1-00 












c 


* 


LU>- 


r— 1 


— 








c 


4s 


o< 


< 


— 




45 




--•* 


* 


CLO 


X 


CO 


o 


— + 




w^» 


+- 


<* 


Cl 


* 


CD 


coO 




X- •4’ 


4* 


K - 


-J 


41 


_l 


* x- 


O 


CL 45 


4© 


vO 


< 


o 


o 


* Ol 


o 


o* 


4« 


00 0 


►. 


X 


*■ 


CD < 


© 


00 o 


4r 


coco— 


o 


CL 


r— 


20 


<\J 


ox 


* 


<r-4X 


► 


w 


— 


<1 — 




—a: 


* 


2U-- 


O 


— — 4 


— 




— 



Vw * 


- - <S> 








— o o 






x=> 


CO — 


O'- * 


- II ^ X- 








-o-c 






— 1 


* CL 


OX # 


OC* Cl 






c 


O c — ' © 






4 — 


* < 


©~ * 


inOcsj < 






o 


Q^ 3 M 






O—i 


— CO 


4 f 


X- — * O 






0 


© 1 1 ^ 






O r*! 


r-i *— i 


— CD 45 


•-cl# 






CO 


roo— o 






©o 


XX 


X * 


- UJ>- CO 






X 


XOOO 






HO 




4 - CL 45 - 


ooCL c _i 






o 


o «o © 






— 1 


200 


— w # 


o- e x- 






o 


om ©ro 






X r -4 


1-0 


02 •» 


*: ►X<CL 






© 


e — CO— 






— r-l 


OOO 


-K 


- ocno< 






<<> 


m 4 ?- 4 - * 






— U_ 


oo 


©00 * 


►in*- uju 






<— - 


* 






COO 


* N — 


000 ^ 45 


'Oh^ 2 - 






* 


* sfCJ-> 






4 f — 


— — 


x< *- 


© ►coCO 






4 © 


41 - * ■«■ — 






41 *- 


COO 


— O * 


CO- Slooi— 






*»» 


— * 4 r CM 






OO 


* r-l CL 


— K * 


—to*: -cl 






— 


— oo* 




> 


+ 4 - 


M -!< 


cm— * 


UJ S'- — < 






CO 


CNJl- *- 




_J 


O 


OOO CD 


* — * 


- v: -xo 






* 


4 t QC — O 




2 


CO 


— CLX 


* O * 


|| - sD ► ► 


o — 


— 


— — * 


— * — K X 




CD 


oO 


+ oo* 


OX * 


o - e — r -4 


o •— 1 


m 


in vocr-c— oror 






r -4 • 


O 1 O 


ia * 


►I-nOCOI— o 


ov 


>* 


> S-XS-XXO-— 




s: 


— CO 


0 — 0 


a:* * 


CCZ. e— < -CM 


<f ► 


► 


► - 0 C 


-c*: ooo* 




o 


x-x 


© 1-4 © 


— o * 


UJU_COUJCCr-i 


r -4 


m 


in so— 


f^- — coo- 




*— i 


CL — 


1 - 0 — CJ 


JO * — 


a. a — s:> 


Ol- 


»— 


1 - 1 - * c— O — 




»- 


OCJ 


— 4 - — 


O © * X 


►< UJ II c ► 


^-2 


2 


2 20 20 r-i *■ ©3 




n 


— 00 * 


4 fr OX 


Oco 4 x 


oa- i- - j— • 


00 


00 


oo OOO to 0 — 3 CO | 


H 


20 


< 045 


— ►-CL 


© — * - 


► Ii 2 X 0 


O ► 


► 


► ► © 


- « 1 — 1 — 


r-l 


o— 


O * — 


CMCLC 


-ix 45 -*- 


cjcooo <x -o 


OO 


o 


OinOfvOco oocDOO 


o 


CL ©CO 


LL— C — 


* <oo 


— O Hi- 4 * 


s 00*— H- 31 — 


o 


o 


0 X 0 i 


O 1 coco 


o 


UJCO* 


2 0 »— i X 


45 OO 


X-O * * 


X<OOOS- ► 


*— oG cO 


©41 ©O 


©04 — I e o 


► 


OX 4 f 


Ox< — 


a:— x 


CL o ** 


►X 2 Xr -4 






r—i — r— ( O 


r-l C tO Q (O f\J 


f — 1 


CL— O 


oo ox- 2 


2 4 f O 


am * * 


— 0-0 -CL 


o ► 


O -o 


► e 


► © 00 * — 4 - 


f — 1 


CD CM 4 


— OUJ— . 


<oo 


OO — - *• * 


X *U)UOG 


ooooo 


o* o^onc- — — 


LL 


# O 


2 exoo 


-J 1 © 


O — 45 * 


X OOO Io 


oo 


©o © 


o* o— 


C— * 


O 


X- 4 s O 


| C 


X — r -4 


— | — * * 


—ro- >- ocri— — 


O e «— * ©^» 


eO ©— 


c -*o* 




003 © 


OO | —45 • 


^ 4 


— * -* 


— >4 1 — ►► — O 


cO 


o 


OXOx 


UVDOO, 


r— 


CL + r-l 


OO 


4 - •— i OC 


— — o— *»00 ►Omh^oO^ 


4 — + 


— or — c. 


— 0 X 0 O 


>t 


— o — 


O* Z 2 f\J 


— 00*1 


1 — oo- «oocm>ooo 


• >o 


vO *— l 


vfISvOOvCDO ©^-l-r 


o 


LLOX 


«— .<< 4 f OOO CO 


— © 0 — 1 * r -4 CO 1 1 — < 5—1 r- 4 -4 




► 4 HflO 


r-lCr-l © 


r~> ©OCVJ*——* 


H 


fOOXh-h-H *— i © 1 -J 


1 , — CM — * 


► -r-l -CL ► ► ► 


GG>G >-0 


OOOrlOH 4 - |- 


► 


UJ— |Q II M << 3 « 


mhhX 


CvJ'— >£•—•« 


> 0*0 — X.O'OvDOD 


tO 


O vt 


O ©G— 


C— *-*Oooa: 


o 


X*— ©OOOO-i 


00 — rH* 



-h* X- — K* ► - — — CO II O II 

<* CL LUC* UJUJ< - rLLIUJlL'X O 

IOOl-H( hh>'WOKKho:>JrJj^h 
t— »— CC C<V 

JaCao^o'ntoyhaKaiKiKOO 
< — — 3U_* j£ 3.LL • ►3:32hUUUUO 



<cc 



< CD 



0 ^ 0 + C+ — OOO 

— vn f — || e II OO 

-j ii «i>-j>-^^a 

— Ji>o r H -J it '-•——i'-' 

<iS)C ir,<us)iL ■& O v 
UCUCuCUilOD 
< < 



3 II r- 4 ,- 4,-4 II II + CL II U- 3 

cl cl ii t-<fccao ooclo n 
x<clookO<o l <: — cl 

OcC<T O-LULUCC c — co *• <r 
UIcC^njJmhhjDCD 
< OX — X-X4-CL 

H(\l r-l 





> 1 - 


o 


H 




o 


o 


o 


a 


o 


o 


o 


a 


CO 


o>t 


OOvTUUU 



68 



oo 


oo 


o 


o 


OOOOOOOOO 


OOCOOOO 


OOOO 


O 


ITNvO 


r— 1 CNJ 




er 


vor- cr(T QH^r >t 


or O' OHNr^ vf 


vc r- <x, o 


(V 


00 s 


OO 


O 


O 


r- r-« r-i — CM <\J (V (N t <\| 




fr fo rn\m 


v* 


OO 


— — 


r— 1 


r— ^ 




— — — — — 1 — 


r-> t—> r—i—i 


r-» 


oo 


OO 


o 


o 


COOOOOOCO 


OOCOOOO 


OC OO 


o 


CC) GO 


CC CD 


CO 


CO 


cc. cci cc cr. oc cr cu cr tjj 


CC- CT CC CC CC CD CC 


cccc era: 


CD 


a: or 


or or 


ar 


or 


ar or or ar ar or or ar cr 


or or or or or a: or 


cr or or cc 


or 


OO 


OO 


O 


O 


OOOOOOOOO 


ooococc 


OOOO 


O 

















<1 


vD 








— 


















to 


o 








CC 








CD 












CM 








© 








O 










< 


»— I 


— 






m 








UJ 










O- 


X Ui 


O 






1 — < 








X 












» - 


o 


— 




UJ 








c 










z: 


* 11 


e 


CM 




- 








«- 










o 


* — 


»— i 






11 








or 










o 


* — 


| 


* 




(NJ 








< 










X 


* U-l 


«— 


»— ■ 




X 


— 






CD 












* -J 


<— ■ 


t— i 




— 


- 


«— 




3 










cr 


• 


CM 


U) 




<M 


4 


— 




C 










UJ 


to ► 


4 


-J 




4 


¥ 


o 




k- 










»- 


UJ o 


4 


X 




4 


4 


o 




•> 










U- 


— o 


or 


w 




CD 


4 


o 




or 










< 


h* H 


< 


¥ 




X 


to 


C\J 




< 








X 




*-H ► 


CD 


•— 




< 


0 r 


V- 




CD 






— - — ■ 


h- 


X 


o ► m 


-J 


M 




— J 


UJ 


o or 




UJ 






or or 


or 


i- 


Oin e 


X 


UJ 




w 


1 — 


e<i 




X 






<< 


< 


or 


-J 1 — CM 


w 


X 




1 


UJ 


C CO 




»- 






CC CD 


UJ 


< 


UJX r-t 


4 


4 




r-i 


X 


CT)— ' 




► 






*—!•—< 




UJ 


> — UJ 


or 


C 




K 


<J 


« — X 




or 






XX 


O 




mx - 


<i 


O 




3 


or 


•r ~ 




< 






w w 


K 


O 


O PvO II 


CD 


• 




4 


< 


4 VO 




CD 






coco 




H- 


ZCM 0 — 


X 


CM 




ar 


Cl 


COO 




cn 




* 


oo 


UJ 




— r-l(Nj ^-i 


4 


4 




< 




<ro 




CL 




>— 


oo 


> 


UJ 


ouji-i uj 


O 


O 




CD — 


h- 


— CC O 




tLI 




i-cm 


OC 


i— i 


> 


o- U.i — X 


o 


O 




JV 


*— 1 


— C-K 




► 




— 4 


* * 


1— 


t— i 


K II - — - • 


— © 


e 




- ►CD 


cm— — ar 




or 


— 


04 


— — 


< 


h- 


3— II uj ►’ 


XCM 


H 




► ocor 


* — * c 




< 


> 


03 


CO CD 


—1 


< 


02 — 00 


► 4 


w 




X oO 


4 cm3cO 


«— 


CC 




J + 


►OC 


U) 


_J 


- — ^'OO 


1— 




irm 




or* xo 


O 


o 


o 


LLO 


CMUJLLI 


or 


UJ 


C2 ►2—1- 


oO 


or 




r r-i 


h- 


<4 XX 


o 


► 


z 


>o 


¥ XX 




ar 


— ZO- — ►(NJ © 


o 




GO UJ UJ 


COO— * 


c- 


ar 


< 


0 


4 OO 


x: 




— <vj- ►UJvOt 


rH— i 


to 


O 


©- 


O 


0+0-“ 


CvJ 


< 




LL*-< 


3 — — 


o 


z 


2 Hoa ©ujw 


o 


o 


m 11 ar 


— xo oo 


O' 


CD 


— 


O — 


4 ^CO 


o 


O 


► CD ►^ ►rv- 1 — 


i 


s 




< 


tt-C ix 


IT, 


Cl 


X 


1 


0 — 0 


X 


t— 1 


— Z'CK— r~» 


li ar 


O 


(NJ 


UJIK 


<4 60- 




— 




toorocvjc 




K- 


— — o — U.h 


-o 


o 


X 


- 04 


coar— «4 a 


NT 


or 


CO 


k-cto 


eOO 


or. 


oo 


>XCvJvOUJ- i 


— to 


© 


— 


II <4 


I<—<H 


f— * 


<f 


UJ 


ZW c 


►^4 * 


UJ 


ox 


►O—i C — 1 II UJO 


r-i 


CM 


— C4 


* CGI CL r-i 


© 


CD 


— 


UjQh 


—a c. 


i 


4- 


O — OU-JCMXUJO — 


w 


4 


< 0:4 


oxorocoro 


< 


►— 


ZK 1 


4 — — 


UL — 


c\jor 


a — z- h — uj 1 


- 4 


4 


4 


CMK- 0_ 




04 OO 1 ox 


«- 


1—1 


O— <— 




<*H 


a c 


e 3— ii uj— or 


►— 


— 


cc<_j a 




« 0 VO eO 4 o 


or 


o 


d_ — r\j 


— >> 


3tOU-0 


r-i<. ► —OO — 


— 


5: 


ariLKi 




wOCimCh 


o 


<r 


o 


XCV4 


CM | 4 


X'X r-J C 


— 4 1 ~ — -n > ii u ujvt or 


— 


<<10 




— O'fr V c-K 


c 


CD 




C Sr 4 


■ CD — 


tom 


O — < 2 4 - hXGK< 


— 


-JO- 1 h- 




»— rv ar o r-i 4 o 


o 


UJ 


o-K or* ccccor^o 


r-<-J rCO — O ►— — 


► CC 


IU 


— X 


►<rto 




CC4<D 


CO 


X 


> 


ar<.oroo<h-c 


©<oo 


e> —4 cO> OtOvO X 


X 




nr uj 





CC 

< 

CD 

X 



(30 CD ©—-Or-* 

4 C — 

O o X ©4 

X.-I I 4c OCM0C 
3 — OX<x 
| hZcif oOCD'*«“ 

a; oc' <r <t cm c uj o — 
<oi— crv cXmm 

CDm <<7 3 CD|- CGCC 
JCO-> v — 

X || (I tti OL -tr H vOvD 

ii arcroro o'- c 

ct«i<7<ia:a:<iLuuj 
< CCjCCCCI/XX CCil-h 
CD 00 tU ZD Ct COLL. — — 
OQ.x<uj<iQra: 
UJ I — I — * 4 |-22 

r-i(NI 



a: r-i l <s co 
<.0— CTCO 

cco»~<o_j 
o I * Q-wo 



0 <cc<ujucdo ©r^or | c-i* >->•* eg- — a ©-© 4 

o zcd^jccs: s: — id— xo h+-o ►• k ►to - (mo o 
2" Mwxwccxa + ciL-M moc« ►►mo. -hdoc 

— ca — a — wwm^CKm + mOC 0DX'0(\'-OLU e4 e 

x ouu^-uji ^ zrzo— <!Z> > © — o ckh^* cm 3 <njo- *-« 3— 3 

i_j ►— ~lj— — > o— 4 • — — t-o—(\j ►ok- ii — x ^ 4- cv— + co- - 

O 3XCMXOtO(/- Jw® a:*— 4. 4* UJO ►r-tvCO ►wVXVi-OII.ClO ► ►- 

z co^ ococ<whz<>mmoocuj er'CMO<c-ou<rooo ► 

— 04 04 4 4 — Z 4 <C 0 ►iLUD^Qll'fOw^fOIIx 

O LLi eX eO- Q. C-K<T , ’* ,M X II D>< ►«—( || — • sD — a o <t © O ► 4 — J ►» — i «-4— 

cr + cm — — — whvc^i-u.ww. hC th ii oxxo—- • 

i- *CD-ow‘-HMi:<Q.uj ii — ii i -hD* uih *-o-i-ii ii — »— i— 

<3 — cQ. II O II DD> m D II — I *— ■ II I- LUO II h<C II II II rsZUJ— <NJUJ<<< 
tcax I oio. can h n n ca — id— — Zk~ x ►Mt-.iH Q :Q:MMK«Krzz 
— X ii ii it O' - l,— 'lu*— *.u»— i>— '►—• culo*— <o£C£h- < < r — >— — ar cr — or or or 
x cj — • — 1 — — — Xtoto:*rujxujuucrarc»- 32 CxcDrrujujftr<<a:ooO 

►UD>2 3 — >3>2»~aa X c 



2 U_ ►— 


u. ► 


OOOC2a.a2U-U.U- 


r-ICM 


<r 




O 




CMC 


O 


*— 1 


oo om 


O 


O 


OCCO 


cn 


0 


O CD 




CD 


CD 



69 



oo 


oococo 


m4* 


4>XC0C Oro 


4-4 


4 444d(U> 


HH 




oo 


oooooo 


CDCC 


cccr cocccccd 


a: a: 


do. or ora or 


OO 


OOGCCO 



cco oooooococoooooooooooooooooooooco 

tnmx- a h atm 4- in 4>r- or O'Cntym 4- in 45 x m O' oh nw 4-in«n xcnc' oh 
mm in m4>4>4>4>4>4>4>4>45xxxxxxf'-xrrncncca: occccrcrcroccoCT'O' 

CCO OOOC C-OOC OOOC'OOC c COGOCCOCOOC ocooo 
cccccD cccoccccccffiircDcQcccrcccccCccccQicCffitcccciritCffioCjCCcDcDtt coaj 
a: or aaaaaraora'ara'aaaofararaararaararaororaciraaraaa: 

OCO COCCGOOC'OOCDOCOOCOCDOCLOOQCOCCCO 



x 


V 




•M 




e 


















4) 


*. 




* 




o 








z 










o 


O 




A 




« 








CD 










CM 


c 




A 




e 








*— i 






PH. 


LU 


pH 


(M 




A 




e 




o 




h- 


o 




X 


a 


111 


rH 




A 




0 




X 




© cO 


UJ 




w 


X 


* ► 


— LU 




A 




0 








JJZ 


X 




H 


X 


II vO 


V- 




A 




0 




o 




ccx 


-J 




O 


c 


a ► © 


II 




A 




0 




1 - 




>>u- 


c 




IL 


a 


C*V<\J 


• H O 




A 




e 








aa 


> 






o 


CC vOrH 


sev in 




A 




0 




-J 




UIUJZ 






z 


c 


C rUJ 


V.X! ► X 




A 




c 




X 




l-h-C 


-J 




c 


x 


<\|- 


vh 4) ► 




A 




e 








ZZm 


c 




a 


G 


• pH || 


o O 




A 




e 




z 




•— 'phOO 


a 




O 




►UJCO 


v- CvJ c 




A 




o 




o 




P— 


o 




o 


00 


c- c 


v* ►h nj 




A 




0 




or 


H 


UJ UJ o 


UJ 




a 


00 


in u o 


h 








o 




LL 


Z 


IIUJ 


i — 




CL 


X 


1— a uj 


*VvO» S UJ 




z 




0 






UJ 


h ►— a 


z 




O CD 


c 


►< 5 : 


VK II » • 




o 




© 




X 


z 


Cl 


H 




LUX 


o 


4) CD CD 


"X pX4> II 


— 


PH 




e 






UJ 


LL LL 






a. oo 




©oo» 


— O 0- c — 


• 


X 


— 


e 




or 


K 


COUJ 


z 




PH 


X 


cvn ► 


00 UJ c<(VS 




C 


m 


© 




LLI 


c 


-J 


o 




xz 


z 


HLLQ 


SCTMUH ►• 


UJ 


-J 


o 


© 




> 


1- 


O CD CD 


PH 




oo 


pH 


hi in 


men* luo ► 


a 


X 


z 


e 




o 


OO 


zzx 


OO 




LU — 


o 


• • X 


pH LL ►• oO 


c 


o 


►H 


e 








XX o 


— 




ari— 


Cl 


ii ► ► 


UJO0*- If || f\lf\J 




-J 


a 


o 




o 


—J 


OC.O 


o 




o 


1 


oro'O 


• >11 4>HHX 


oo 


c 


CL 


© 




UJ 


c 


CDCC 


UJ 




OOZ 


4) 


c in o 


II CrHKCUJ ► 


a: 


o 


o 


e 




z 


z 


-J 


a 




z x 


H 


cCh-nj 


Jero pi* \ 


UJ 




X 


c 




X 


a: 


oror< 


CL 




CLL 




-J •- H 


ujoovoa n ► 


X 


to 


O0 


© 




^5 


UJ 


Uu LL' Z 






a © 


LL 


OIL! 


CLUX c«JX4> 


u. 


ph 




© 




O0 


h- 


Saar 


LU 




CZor 


O 


m 0 » 


K» N<2 c 


2 


X 


O 


• 






X 


DO. UJ 


—J 




DDL 




►cm ll 


* Z ►pH- CM 


c 


X 


c 


0 




«- 


LU 


__J DD 1 — 


CD 




a— oo 


00 


Oh or 


*■— in uj ►• p-i 


a 




0 


© 




— . 




X 


X 




a oox 


z 


COLLI C 


o m • in ►ui 


c 


a 


pH 


o 




X 


z 


ZZ LU 


oO 




CD PH 


c 


X- CD 


inujx ii voo- 


a 


o 




© 




w 


— c 


OO 1 


OO 




X OIL’ 


UJ 


► no 


xx ►nj|- m n 




IL 


z 


© 




h- 


X 


mmZOJ 




00 LUX 


z 


o arc 


►1 — vO^H .-XH 


_j 




c 


e 




O 


►oo 


000000 C00O 




an 




eC x 


00 ©X4) ►(— 


c 


-J 


X 


e 




LL 


l-U! 


a PH- 


xz 




Z0 l 


X 


roc c 


eVMQ. o4> 


►— 




X 


• 




w 


Oar 


IL u GIL 


— 




O X 


CD 


hC» 


mvHCN o • 


X 


c 




© 




-J 


LL i— i 


L-lLUJ02h 




P- LU CD 




UJ ► 


HVUJU-HfM h 


ph 


V 


a 


« 




C 


►X 


Luaa ■ 


C-J 






UJ 


- • o 


ID’X- - LL‘HO 


Z 


V 


UJ 


© 




a 


DO 


ZQ.aii'orx 




occo 


z 


II ►CO 


• v ll ►- a jin 


ph 


V 


X 


© 




o 


XLU 


C ZOoo 




ZDOJ 


o 


a ox — 


II SHlf || • X 




V 


C 


9 




UJ 


POT 


aUJlUCCLL' 




xox 


o 


cm pv 


2 VD-rGll ► 


UJ 


V 


UJ 


e 




H 


—1 


cjjzcror 




ILCOO 




CDX4) * 


<vci-. o o 


I 


V 


or 


V 


vO 


Z 


XI- 


CL CD CD 1 


a 




— 


00 


X ► o'COVO ► ►! o 


X 


V 


o 


0 


pH 


•— i 


p-O 


XXUJCQUJ 




OJZ 


PH 


4>CM c 


i xx4'anorc\i 




V 




© 


O 




u. 


ILODIDI 




zca 




• CrHC\j 


|- - OfO» pH' 


V- 


V 


o 


0 


G 


UJ 


4) 


CCCI-oOl— 




czx 


z 


►CMUJfH 


i . ►rvx puj 


*■ p 


V 


X 


« 


O 


1- 


ho: 








au_ 


c 


Oh. LLJ 


i ► ►Oh pQ. 


sOO 


V 


or 


• 




X 


OUJ 


Z 1 | 1 


l 




O0UJ 


PH 


HLU |l - 


OOhLUvOCM II 


cm 


V 




© 


UJ 


CL 


Ch 


o 






UJh-LU 


X 


X» all HCMX- cX < (\JX — 


• • 


— * p 


© 


z 


2 


OUJ 


PH 






ZXQC 


c 


► ll caxx ► ii cv ►c. 


PH -IT 




pH P- 


© 


•— 1 


O 


z 


H 






— UJ 


X 


Oa ccc 


► ►COHOCCLJCO'O 


ao 


• 


H- 


IUO 


JC 


a 




OOIL' 


►— K 


_J 


icaiaxxihtiis. 


► 10- 


X 


crx 


• 


O 


00 


jo: 


— »- 




xz 


XUjoo 


cc 


HCDIM 


hhhQ.. hC l| h • 




► H 


© 


o 


oo 


LU C C 


otldo 




aro 


CXX 


c> 


-CLh* 


< ii -'_iph— n 




© 


or 


CLI- 


CUQ- 


uxxa 




cz 


ax z 


XU) 


1 - 


■XXXOXX 


CO- W 


X 


— l- 


« 


CD 


a 


< 


00 




z 


CD 


X 


c» - occc- • c- 


:cu 


CCLUC 


c 


D 


X 


oo 


LU 




UJ 


X 


LU 


21 ► ► HS'Z Z ► *2 ► 


m 2 »-2 xx 2 a 


0 


oo 


a 


X 


c 




a 


oo 


z 



or o ox ar a or o O or o 3 or — or —a oc o 



Chh -CCChhON GC'Crt^Ct-Z e 

LLXI ^LLLLU. U.2il02U.^H1 o 

< CCO CCO CcC 

■4 O' 

CO O' 



in 

00 



HCjm H 

oo o 

H H iH 



O m c»-*pv 

o o a O' 

h O a 0 s OOOOOOOOOOOOOOOOOOOOOOOOOOOO 



70 



ooooooo 
<vrn<t-ir>vO^ oo 
a a O’ o o a a 

rH r— *1 ^-4 r— 4 r— 4 r— 4 

ooooooo 
cr cl cd cd a cc cc 

CL CCcCCCCC DCCC 

OOOOOOO 



LUO. e 

laiQ .3 « 

O < LU © 

♦-a: Z c 

XO •* O 



ooooooocoocooooooocoooooo 

a o rococo m- in soc-cc o-OHfvcovd-irivOh cr o Or^njco 

O'OOOOOOOOCCHHf-lHHHr-t^rV^^'Mrj 

*-4<vcmcmcvcmo<c\ cmcmoocmcmcmcmcvCMcmcmcocmcmcvcmcm 

coccccooooooooccc cccococo 
cccccLaccrcccLcCcDCDcr cccraDcr cdcdcdcd cd core co clod 
ccccccccccacccccccxccci'cccccccccccccrcLcccc-CCCCcC 
OOOOCOOCC CC-CC OOCOCCCOCOOO 



3LL> •* 


c 


>0 


Z 


OX 


• 


H 


O 


- #o< 


0 


O 


t-rf 




ft 


O 


CO 


JJOCJ 


e 


O 


*-H 




c 




O 


zo 


• 


LU 


LU 


cc <*-^o 


0 


z 


CC 


ox ©< 


ft 


*—> 


CL 


LLLL>2I 


c 


k- 






ft 


3 


UJ 




ft 


O 


-J 




ft 


a: 


CO 




ft 


CD 


3 




ft 


3 


O 




e 


OO 


O 



oooooooooc oo 

LT\ vT C- co a C l/> 

CMC'JCocMCvcvcoror ri roroc^ 
CM fV CO CO CM CO CO CM CO CM CM CO 
OOCOCoOOOCCO 
cc cd cd cd cd cd aj cr. cr; cn cc cr 
ccccccocgcccccccolcccccl 
COCGCCCCOOOC 

* 

o 





0 


ft 
























X 






0 


ft 






















3 


a: — k- 








ft 






















*> 


mcl 






CO 


0 






















h- 


* —a 




r-4 


1 


ft 






















CL 


km 




co 


I s - 


0 






















< 


0 cco 






CO 


e 






















O 


0 0 * 


1 


LU 


CO 


© 






















•* 


ft— coc 


NO 


LU 




e 






















CO 


CM— OX 


O 


nr 


CCO 


© 






















r-4 


— —'CL 


O 


O 


JZ 


e 






















k- 


1 04 — 


K 


UJ 


<< 


0 






















CL 


0 z 


•— 


O 


or 


0 






















< 


— c— *— 


CO 




Ch 


0 






















O 


— CMOCO 


0 


O 


LLI 1 — 1 


c 






















•- 


- — "Oor 


0 


k- 


L— r-4 


0 




















— 


r-4 


CMO c< 


0 




z 1 


0 






















O 


* C (MO 


* 


Cl 


—0 


0 




















-«■ 


CM 


* 0^* 


KJ 


3 


0 


« 








— 




— 








O 


r-4 


NJC’I — 


* 




LL «—* 


0 




















1 


> 


W w^- — 


O *~ 


CO 


O O 


e 






•— 


O 


O 


O 


— * 




— 


<r 


•* 


** MO 


O — 


0 


Q. 


c 






■ — 


1 


1 


1 


O 


0 


O 


«— 


r4 


— * * X 


ft — 


< 


ZQ- 


o 






O 


<1 


< 


<f 


1 


1 


1 


»- 


O 


CMO— or 


CM M 


t— < 


O 


e 






1 


w 


w 


w 


< 


< 


< 


0 


O 


* CCM* 


1 1 


X 


►— < V 


e 






C 


k- 


k- 


K 


w 


w 




LL 


h- 


* ft * 0 


O 


0 


KCO 


0 






w 


0 


O 


0 


1- 


k- 


k- 


4- 




Ccm* Q 


(MOO 


z 


<nO 


0 






k- 


LL 


IL 


LL 


0 


0 


O 




r-4 


X-'O 0 


* k- 0 


> 


-JO 


0 






O 


4 


+ 


+ 


UL 


UL 


LL 


O 


CL 


or* x co 


* a»-4 


0 


3 r-4 


0 






LL 


— 


<— 




+ 


4- 


4- 


4 


O 


— or^ 


Nx.Kt—' 


0 


O 


ft 




k- 


4 


O 


O 


O 


— 






< 


O 


| (Ow-v. 


wU^ 


CL LU O •* 


ft 




O 


— 


4 


■f 


+ 


O 


O 


O 




K 


CO* 1 0 


4 1 O 


UJCZ 


0 




LL 


O 


<t 


< 


< 


4- 


4- 


4- 


>- 


•> 


-'*00 


OfM-' 


CO CO OO 


0 




*. 


4 


w 


W ' 


w 


<1 


< 


< 


O 


O 


HOG o 


COM 


UJ 


O 


e 




O 


< 


k- 


1- 


k- 


w 


w 


w 


LL 


•> 


CC X rCO 


ftOO 


K 




0 




*> 


w 


0 


0 


O 


K 


k- 


H 


— - 


r-4 


OCCrHw 


— r-t k-C 


cu-'k-o 


© 


'> 


CC 


k- 


LL 


LL 


U_ 


O 


O 


O 


* 


< 




0 


DC O <1 — 1 


0 


*> 


•> 


0 


w 


«w 


w 


U- 


UL 


LL 


H 


X 


Owl- 1 




ezi\ 


0 


k- 


< 


LL 


* 


* 


* 


■W 


--r 


w 


1 


— CL 


* or 




LUUj^^ 


0 


O 


•> 


w 


r-4 


»-4 


r — 1 


* 


4< 


* CD O 


4^-J 


0 — 0 — 


ft — O 4 


haxa' 


• 


LL 


> 


CD* CD I 


CD 1 1 


CD I CD r— ICC CC 




* IO 


r< 


X — CO- 


CMN.G- 


ZujCO 


e 


•> 




* -4* O* O* O* 


1 * 


1 * 


1 r— 1 CM 


M r 


OT— O — 


w 1 | N 


*-LL(X> 


• 


3 


3 


O 1 O'COCVO'COCOQCO 


1 


1 O 


* CM* isll 


h- — — < 



coc^caoroh or-oo'orc 
c-r~ O'O'k-coroo c-coocMODr-C'io 

OvT c.'cccc'Oir\vOcooc<j^<v} r-^coc- 
»tO'OC\iir>Nl"r-ivOf-tOO<COO''Oh-CO 
c\ic-»-(cCr-ic\JOC'Coc\jr-Hircv | vj oocv 
COh-sOO'O'^-UPC-CCCOODf^NOCOrHlC 

irscO'Ovf>rc<t'f^-fvjf-( 00 oO'OrMccO 
— oincco'O'cor-'OoNf cosD^i^cocc 
_j ccorc.r^^in 

X C-r- r-l'Or-ilO nJ- a-’ CM O'er) ITMTNOO CO 

+ oDo -^ocMo>fr coO'cccvr~-f^ ^ r- 

3 stcoic\' 0'0 0 cvj't*-»f s ^cccch'*^iriv} 

X OCOh-CMlAr-CM^OOO'CCr-^OCvjO' 

— - OvOcor-4r-4inococoh-oioocovO © 

* _Jh-r^c\jr--icnr^r^c\jO'<to-ocr'^o-»- 

OX'J it cMcoc\i-o-h--vOcnr^c cc go io>- 3 

O I OccC- ©CO # h- (C c CM csj- Of^— 'DC 

in 3 >t —« 4 M 4 ro 4 cn 4 (M 4 «- 44 >l-* 3 
eX e o ©>- ©>- ©>- ft>- ©> o> rtThD 
II II II II II II II II II II II 11 |l II II II II IIUJZ 
<CD0>0>U>0>'0>0>-U>-0>Cr:iU 



O* — * nrc— x 

, o* — — c?k-— a. 

r eNNwi/iQ-O J 
} CO w--* 0<G< 

► — V k-* * O 0 1 

3 Z — OTC — | CO r-4 

: *-«cooxvi“fvi'v.< 

C 00 * LOCL* OOX 
< a:* o—* oca 
. < 0 — 1 ck u 

J ox-oi- cc< 
3 wCtr^CDtToO^— ' 
3 — e ~ 0 * * 

I I I — •-'O* CO 
«. NO II — ^O — O 
I II^NI LL 

3^ k-<or oil 

' 4 NaXO- li l — 
'ZOCcloo'-mcu- 
3 II CWJDMC<Z 
;ZhC<^U. OcO 
< < CD < 



ooooooo 



71 



oo oooooooooooooo 

\OC* cr crOr-*<vrovt lOvor^coO'OrH 

rcro vj vt in lo 10 tr> ir trio ir> it LT ^o vO 

COCO Cv, (N : rv CO CO CO (VI CO Cyifv TV CO (V 

oo OOOOCOOOC .OOOOO 

cnco ODCCCCCCCDCDCLCCcf.CCCCcCCCia) 

ex ex' excxcxcxcxexexo'cxcxexcxexex 

OO OCOCCCGCOCCCCO 



2 

£X 

3 

t-O 

u-'zr 

CXLU 



2 


* 


ft 


O * 


o 


X * 


ft 


ex — 


K 




Cl 


* 0 


< 


►— 


O 


0 a 


ft 


0 < 


CO 


C — O 


»— < 


CO— - 1 


1 - 


— ' NO — 


a 


1 00 — 


< 


CO -* 


0 


— eh- CO 
— <v— | 


•> 


— V.O 0 0 


» — ! 


NOD O 


O 


* Oo c 


CO 


* eO *— 1 


H 


iv- w. 


> 


wwPv.' * 


► 


* ft * a 


*— < 


— * 0 *- 


O 


— coco s 


O 


— * 0 © 0 


H 


— * r CM O 


•* 


rv! Oco | 0 


«-» 


*r X — — r-< 


a 


— ex* rvi — 


0 


CO * 


0 


* 1 CO*. + 


H- 


* CO* K 0 


•> 


O W* w — 


O 


X hD+-M 


ft 


ex cXXCO< 


*~t 


— OXCKH 


< 


1 on w ft a u_i 


X 


O C^H<CD 


— Q 


O * v. — 0 1 


*A_J 


e O — — 1 ^ 


►< 


*-4 X VN< 


fvl - 


— CX— — .OK 


1 O 


* * CO C OLL 


CK 


— 0 — 0 * h-t— CG 


►a 


hCNO# CX — w 


x< 


a ea eN<Z 4 i 


1 u 


OMCW-Omn 


<r *■ 


0— 0— > 1 ono 


— — O 


KM-Z-10CO 


CJCCI 


\.W>--(\,Q * 


-'ft or 


OwOon* C 1 CO 


rc_iv. 


CKOex* 1 Cl 


I-<UL 


eCT 0 <K_— — O 


ZOJJ 


hGhCIZ-O 


oner O 


wiow wex^-oi- 


O 


I C | wi/OO* 


Zh< 


1 C eM 


O*— v. 


II It |l CO — CvJ 


►-UZ 


II — *V. 11 ; 




K II 0 < 



o_jx + <a<jfsjcx cr»D 
zarZi-c^cPN ob-v-G 
O 5: C' II LLI C LL C LC o C' Z! U t 2 

u.»-uzcK-iXhCu-ioauj 
< < 



72 



REFERENCES 



1. Szebehely, V. G., Theory of Orbits , p. v, Academic 

Press, 1967. 

2. Lagerstrom, P. and Kovorkian, J. , "Earth-to-Moon 

Trajectories in the Restricted Three-Body Problem,” 
Journal do Mecanique , V. 2, p. 189-218, June 1963. 

3. Lagerstrom, P. and Kevorkian, J., M Earth-to-Moon 

Trajectories with Minimal Energy," Journal de 
Macanique , V. 2, Dec. 1963 . 

*4 . Van Dyke, M. , Perturbation Methods in Fluid Mechanics , 
Academic Press , 1964 . 

5 . Arne s , W . , Non Linear Ordinary Differential Equations 

in Transport Processes, p. 172-179 * Academic Press, 
1968 . 

6. Shi, Y. , and Eckstein, M. , "Uniformly Valid Asymptotic 

Solutions of Non-Planar Earth-to-Moon Trajectories 
in the Restricted Four-Body Problem," The Astronomical 
Journal , V. 72, p. 685-701, Aug. 1967- 

7. Lagerstrom, P. , and Kevorkian, J., Numerical Aspects 

of Uniformly Valid Asymptotic Approximations for a 
Class of Trajectories in the Restricted Three-Body 
Problem," Progress in Astronautics and Aeronautics , 

V. 14, Academic Press, 1964. 

8. Kevorkian, J. and Brachet, G., "Numerical Analysis of 

the Asymptotic Solution for Earth-to-Moon Trajec- 
tories," AIAA J., V. 7 ? p. 885 - 889 , May 1969. 

9. Kevorkian, J. and Lagerstrom, P, , "Non-Planar Earth- 

to-Moon Trajectories in the Restricted Three-Body 
Problem," AIAA J., V. 4, p. 149-152, Jan 1966. 



73 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Prof. D. J. Collins, Code 57Co 1 

Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

4. Chairman, Department of Aeronautics 1 

Naval Postgraduate School 

Monterey, California 93940 

5. LT(jg) Vernon Curtis Gordon, USN 1 

Route 7, Box 230 

Athens, Alabama 35611 



74 



UNCLASSIFIED 

Seouri t y Classifc • at ion 



DOCUMENT CONTROL DATA - R & D 

(Security classification of title, body of abstract and indexing annotation must be entered when the overall report Is classified) 



i. Originating activity (Corporate author) 

Naval Postgraduate School 
Monterey, California 939^0 



2a. REPORT SECURITY CLASSIFICATION 

Unclassified 



2b. GROUP 



3 REPORT TITLE 



SINGULAR PERTURBATION PROBLEMS IN EARTH-MOON SPACE 



4. DESCRIPTIVE NOTES (Type of report and, inclusive dates) 

Master’s Thesis; June 1970 



5- AuTHOR(S) (First name, middle initial, le»f name) 

Vernon Curtis Gordon, Lieutenant (junior grade). United States Navy 



6 REPORT O A TE 

June 1970 


7a. TOTAL NO. OF PAGES 

75 


76. NO. OF REFS 

9 


8a. CONTRACT OR GRANT NO. 
6. PROJECT NO. 

c. 

d. 


9a. ORIGINATOR'S REPORT NUM BER(S) 


9b. OTHER REPORT NO(S) (Any other numbers that may be aa signed 
this report) 


IO. DISTRIBUTION STATEMENT 

This document has been approved for public release and sale; 
its distribution is unlimited. 


11. SUPPLEMENTARY notes 


12. SPONSORING MILITARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 939^0 



13. ABSTRACT 



A program was written for the IBM 360 system of the 
Naval Postgraduate School Computer Facility to apply the 
methods of singular perturbation theory to earth-moon 
traj ectories . 

Verification of results against previous work was 
carried out and it was found that agreement could be at- 
tained only for energy parameter values of 0.707 or less. 
No solution for higher values could be found. 

The analysis of three dimensional orbits was then 
conducted within this restricted range to show the merit 
of singular perturbation theory as an initial design tool. 



DD , F r.!473 (PAGE 1) 

S/N 01 01 -807-681 1 



UNCLASSIFIED 



75 



Security Classification 



UNCLASSIFIED 

Security Classification 



CEY WORDS 



PERTURBATIONS 

EARTH-MOON SYSTEM 

MATCHED ASYMPTOTIC EXPANSIONS 

TRAJECTORIES 



ROLE W T 



ROLE W T 



1473 < BACK > 



01 * 807-6821 



UNCLASSIFIED 

Security Classification 



Thesi s 

G576 

C.l 



j; on n p / 
6 1 1 ^ 



Gordon 

Drobl* 9U,a - r pertur batior 
Problems in earth 
space. earth -moon 




Thesi s 
G576 
r .1 



Gordon 

Singular perturbation 
problems in earth-moon 
space . 



thesG576 

Singular perturbation problems in earth 




3 2768 002 13136 9 

DUDLEY KNOX LIBRARY 



